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We study the response of the QCD vacuum to strong magnetic fields, using a 
potential model for the quark-antiquark interaction. We find that production 
of spin-polarized uu pairs is energetically favorable for fields B > Bcru ~ 
10 GeV^. We contrast the resulting uu condensate with the quark condensate 
which is present at zero magnetic field, and we estimate the corresponding 
magnetization as a function of B. 
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1 Introduction 



Strong magnetic fields are interesting from several perspectives. From a 
theoretical point of view, an external magnetic field allows one to probe the 
vacuum structure and correlation functions of a quantum field theory. Strong 
magnetic fields are also of interest in astrophysics. There may be neutron 
stars with fields of up to 10^^ - 10^^ Gauss while it has been suggested 
that much larger fields existed in the early universe 

With increasing magnetic field, the first place one might expect something 
interesting to happen is at the scale set by the electron mass, 

5 = m^/v/^ = 4.4 X 10^^ Gauss . (1) 

At this scale an electron's Landau energy equals its rest energy. Magnetic 
fields of this strength have a significant effect on atomic and molecular 
physics, as reviewed in 0. However the structure of the QED vacuum does 
not change dramatically in fields of this magnitude. Corrections to the en- 
ergy of a free electron in the lowest Landau level are small, proportional to 
the electron's anomalous magnetic moment, so it is not energetically favor- 
able to produce e'^e~ pairs. The binding energy of positronium is small and 
does not change this conclusion, so the QED vacuum is stable.^ 

The next place one might expect something interesting to happen is at 
the QCD scale, 

B = Alco/Va ^ 10^^ Gauss . (2) 

This regime will be the main focus of this paper. We will argue that the per- 
turbative QCD vacuum becomes unstable to formation of a quark-antiquark 
condensate. The basic physics is easy to understand: the strong magnetic 
field restricts quarks and antiquarks to move in one dimension, and the 
strongly attractive QCD potential then leads to formation of a bound state 
with negative energy. We will argue that at sufficiently large magnetic fields 
the effective coupling becomes weak and perturbative QCD can be used. 

The next interesting regime starts at the electroweak scale, 

B = ml, /v^ ^ 10^^ Gauss . (3) 

^This is discussed in more detail in Exponentially large magnetic fields, in contrast, 
have been shown to catalyze chiral symmetry breaking in QED ^ . 
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A field of tliis magnitude lias been argued to drive electroweak symmetry 
restoration [Q. Finally, a grand unified theory with magnetic monopoles of 
mass Mmon could provide an absolute upper bound on a possible magnetic 
field. Extrapolation of semiclassical calculations p suggests that monopole 
pair production would become copious and short out any existing magnetic 
field when 



B - v^ML, ~ 10^^ [w^) ' Ga^ss • (4) 

Although the approximations underlying the semiclassical calculation break 
down before such fields are reached, this is probably an overestimate of the 
maximum possible field. Arguments similar to those leading to Eq. @ in- 
dicate that there should be a local restoration of the GUT symmetry when 
B ~ a^/^M^Q^. (The resulting regions of symmetric vacuum can be viewed 
as condensed monopole and antimonopole cores.) With the unbroken sym- 
metry enlarged to a simple non-Abelian group, magnetic flux is not conserved 
and a coherent long-range magnetic field can no longer be sustained. 

There are several approaches one could adopt for studying QCD in a 
strong magnetic field. At the hadronic level, one can study the effect of a 
magnetic field on hadron spectra and the nuclear equation of state Q, 
based on the large anomalous magnetic moments of hadrons. Alternatively, 
one can take a diagrammatic approach, and compute the quark condensate 
in a large magnetic field by resumming diagrams . The effect of a magnetic 



field has also been investigated in the Nambu-Jona-Lasinio model [p!0| , [11], [12 
and in an instanton-inspired model for chiral symmetry breaking [11] . In 
contrast, we study the problem using the quark model. The advantages of 
this approach are simplicity and a clear physical picture of the QCD vacuum. 

Throughout this paper we set h = c = 1. The conversion factor is 
1 GeV^/(/?.c)^/^ = 1.44 X 10^^ Gauss. In the introduction we have used 
Gaussian units, but in the remainder of this paper we will exclusively use 

Heaviside-LorentZ units: ^Gaussian = ViTr^Heaviside-Lorentz and gcaussian = 

9Heaviside-Lorentz/\/47r. Thus, for example, in the remainder of this paper the 
charge of an up quark is g = | ■ A/iTra. 

An outline of this paper is as follows. In section 2 we study the behavior 
of a qq pair in a strong magnetic field, with the help of a potential model 
for the quark- ant iquark interaction. In section 3 we estimate the strength 
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of the qq condensate in magnetic fields somewliat above the QCD scale. 
In section 4 we study the condensate in the regime of large fields, where 
perturbative QCD is applicable. Section 5 contains a summary and some 
concluding comments. For completeness, in the appendix we compute the 
response of the QCD vacuum to weak magnetic fields by performing a pion 
loop calculation. 



2 Mesons in a strong magnetic field 

In a strong magnetic field quarks follow Landau orbits in the directions trans- 
verse to the magnetic field. These have a characteristic radius R = 1/ y/qB, 
so that for B > 1 GeV^ the quarks can be localized in the two transverse 
directions to distances shorter than the QCD scale. 

Moreover, there is no energy cost associated with this localization. Intu- 
itively, this is because the quark kinetic energy is canceled when the magnetic 
moment of the quark lines up with the magnetic field. More precisely, the 
energy levels for a Dirac particle in a background magnetic field are 



Here n = 0, 1, 2, . . . labels the Landau levels, a = ±1 specifies the spin 
orientation, and pz is the momentum in the z direction. Thus in the lowest 
Landau level, with an appropriate spin orientation, the quark behaves just 
like a relativistic particle in 1 + 1 dimensions. 

One might expect that this localization enhances the attraction between a 
color-singlet quark and antiquark to the point where the energy of a qq state 
becomes negative. This would signal an instability with respect to formation 
of a spin-polarized qq condensate. 

To address this issue, we wish to estimate the energy of a qq state in a 
strong magnetic field. We do this by adopting a potential model for the qq 
interaction That is, we will take the Hamiltonian for a qq state to be 
given by the quasi-relativistic expression 






(6) 
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A wide variety of potentials have been discussed in the hterature; we will use 
the Cornell potential [0 



V{r) = Ar-- + C. (7) 

r 



We will be focusing on uu or dd states, and so use parameters ||15 



A = 0.203 GeV^ 

K = 0.437 (8) 

m = 0.150 GeV 

C = -0.599 GeV 

chosen to fit the the spectrum of light mesons. 

In a strong magnetic field this three-dimensional model should reduce to 
an effective one- dimensional problem. However, at distances shorter than the 
magnetic length R the problem again becomes three-dimensional. We can 
take this into account by cutting off our one- dimensional potential at short 
distances. Thus we study the one-dimensional problem 

H = 2y/pl + m'^ + V{z) (9) 



Viz) 



AR-^ + C \z\<R 
R 

Az h C \z\> R. 

z 



By considering a Gaussian trial wavefunction, one can easily see that as 
^ the spectrum of this Hamiltonian is unbounded from below. 

To estimate the energy levels of the Hamiltonian (|^) we use a WKB 
approximation ||16|. The classical turning points are aX z = ±L, where 



V{L) = E — 2m. The WKB quantization condition §Pzdz = 2'n{n + |) 
becomes 



^ dz^{V{z) - Ef - 4m2 = TT + n = 0,l,2. 



(10) 
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Figure 1: Ground state energy as a function of qB. E is in units of GeV, qB 
is in units of GeV^. 
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Figure 3: L/R as a. function of qB (measured in units of GeV^). 



The resulting ground state energy is shown in Fig. 1. Note that the energy 
is negative for qB > 2 GeV^. We expect the WKB approximation to give a 
reasonable estimate for the ground state energy in this regime. 

The semiclassical turning point L is shown in Fig. 2. Our reduction to 
one dimension only makes sense if L is large compared to R. As can be seen 
in Fig. 3 this condition is reasonably well satisfied in the regime of interest. 
Note that L decreases as the magnetic field gets bigger. This means that at 
sufficiently large magnetic fields the qq bound state is driven into a short- 
distance regime where perturbative QCD can be applied. This regime is 
discussed in more detail in section 4. 



3 Meson condensation 



When the magnetic field is sufficiently large, qq pairs will start to condense. 
We first consider a single flavor, and discuss condensation of uu pairs. Con- 
densation occurs when qB > 2 GeV^, which for a uu composite means 

2 GeV^ 

B > i?„it ^ , ^ 10 GeV^ . (11) 

The quark magnetic moments line up with the magnetic field, so B is in- 
creased by the formation of the condensate. This would seem to make the 
vacuum unstable, but eventually the uu pairs will start to interact, and this 
effect presumably stabilizes the system. 

Because the uu pairs are color singlets, they will not interact strongly until 
their wavefunctions begin to overlap. Hence, for B > B^rit pair production 
should proceed unimpeded until the density of uu pairs reaches a value of 
roughly 

Once this density is attained, QCD interactions between pairs will tend to 
suppress further growth in the condensate. We will use Eq. (|12D as our es- 
timate for p, although the actual value that emerges from the interplay of 
magnetic and QCD effects will presumably have a somewhat more compli- 
cated dependence on B. 
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Treating the quark and antiquark as elementary Dirac fermions, the mag- 
netic moment of a pair is fi = q/m and the magnetization is 



To evaluate this, note that Fig. 2 shows that for B > Bcru the length L is 
slowly varying, with L ^ 2 GeV^^. The question of what value to use for the 
mass is a bit more subtle. At zero magnetic field one uses constituent quark 
masses m ^ 300 MeV to estimate magnetic moments, although for extremely 
large magnetic fields, where R is very small, a current quark mass may be 
more appropriate. Using the w-quark charge in Eq. (^3|) gives 

/300MeV\ /2 GeV"^\ „ 
M = 0.022 (14) 

Hence, we expect that in the regime B > Bcru the magnetization will be 
small compared to B, so that we are justified in ignoring the back-reaction 
of the magnetization on the strength of the condensate. 

We now consider the effects of the other quark flavors. If there were 
no interaction between quarks of different flavors, extension of the above 
analysis to d quarks would predict that a dd condensate forms at a critical 
field which is twice as large as for uu pairs, and with a magnetization that 
is one quarter as large for a given B. However, the different condensates 
will interact with each other, so that an increase in the condensate of one 
flavor will tend to cause a compensating decrease in the other condensates. 
Because of the relatively large u quark electric charge, a uu condensate is 
energetically favored over dd or ss, while the heavier quarks are suppressed 
by their mass. Hence we expect the condensate to be dominated by uu pairs. 



4 Condensation in the perturbative regime 

When the magnetic field is far above the QCD scale the qq composite is 
driven into a short-distance regime where perturbative QCD can be applied. 
In this section we discuss the magnetization in this perturbative regime. 
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At short distances we should replace the Cornell potential of Eq. (|^ with 
the potential from one-gluon exchange, 

Vir) = -- ^ = ^ fl5) 

^ ' 3 r ~ rlog(l/Ar) ^ ' 

Stt Stt 

A 



3bo llNc - 2Nf 

The WKB quantization condition for this potential is still given by Eq. (10) 



If one makes the approximation of neglecting the quark mass, the WKB 
integral can be evaluated analytically to obtain a relation between the turning 
point L and the radius R = 1/ \/qB. 



1 A 1 TT 

+ Aloglog— = - — + A log log — - - ■ (16) 



log(l/LA) ° °LA log(l/i?A) ^ ^ RA 2 

For R —>■ this, together with Eq. (|13|), gives a magnetization 



M=^^^n^ . (17) 

mirL mrr \ A^ / 

To evaluate the exponent in the last factor we set Nc = 3 and take Nf to 
be the number of quark flavors that are lighter than the mass scale set by 
B. The dependence on Nf is actually rather weak, with any value between 
2 and 6 yielding an exponent of about 0.1. 

Of course these calculations are only valid if the turning point L is small 
enough to trust the potential of Eq. (|1^). This requires magnetic fields that 
are far larger than those we have considered so far. For example, L < ^ 
for a uu composite requires B > 5 x 10^^ A^; although several orders of 
magnitude above the electroweak scale, this is still far below the upper limit 
ofEq. (§. 

We expect that our approximations give the right qualitative behavior of 
the magnetization at strong fields. However an accurate quantitative calcu- 
lation calls for more sophisticated techniques than we have employed here. 
For one thing, the WKB estimate of Eq. ([TBD gets worse as the magnetic field 
increases, since the condition for the vahdity of WKB, \dp~^/dz\ ^ 1, is vi- 
olated near z = R. (The left hand side grows logarithmically as -R — > 0.) A 
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more serious concern is that as i? — we should really treat the qq composite 
in a fully relativistic manner, e.g., by solving a Bethe-Salpeter equation. 

Furthermore, at sufficiently high fields the nonlinearities become impor- 
tant enough that we must take into account the back-reaction of the magneti- 
zation on the condensate. These nonlinearities can arise from several sources. 
First, there is the explicit nonlinearity in Eq. ([T7|), which shows that M/B 
includes a factor that grows as a small power of the magnetic field. Next, we 
expect the effective mass of the quarks to decrease as B grows, also leading 
to an increase in M/B. A third possible source is the corrections to our 
estimate Eq. (JT^) for the density of condensed pairs; these should also give 
an increase in M/ B at stronger fields. 



5 Conclusions 

In this paper we have used a quark model approach to study the behavior 
of QCD in the presence of a strong magnetic field. In the presence of such a 
field the quarks can be localized in the two transverse directions with no cost 
in energy. This enhances the quark-antiquark attraction to such an extent 
that the binding energy can compensate for the mass, thus making uu pair 
production energetically favorable, if B is greater than a critical value of 
about 10 GeV2. 

In the language of field theory, this pair production corresponds to the 
formation of a chiral symmetry breaking uu condensate. Of course, even in 
the absence of a magnetic field, nonperturbative QCD dynamics produce a 
nonzero quark condensate that breaks chiral symmetry. However, the zero- 
field and the high field condensates differ in some significant aspects. 

At zero field the condensate is Lorentz invariant. In particular, (qa'^'^q) = 
0. By contrast, the quark pairs produced by a critical magnetic field are 
polarized along the direction of the magnetic field. For a field directed along 
the 2;-direction, this corresponds to a condensate with (qq) ^ {q<y^^q)-f\ 

The fiavor properties of the two condensates are also quite different. The 
'^The possibility of a {qa^^^q) condensate was raised in 
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zero temperature, zero field condensate is, to a good approximation, flavor 
SU(3) symmetric, with {uu) ~ {dd) (ss). This is not the case in the 
presence of a super-critical magnetic field, since the production of uu pairs 
is energetically favored over that of dd and ss pairs. 

Finally, the zero field and high field condensates differ in magnitude. The 
former is of the order of Aq^j-, ~ (.25 GeV)^. This should be compared with 
our estimate, Eq. ([12D, for the density of quark pairs. This density increases 
faster than linearly with B, but even at the critical field for uu production 
we have p ~ (0.7GeV)l 

We would like to understand the transition between the zero field and high 
field regimes. As B is increased from zero, its initial effect is to gradually 
polarize the QCD chiral condensate. For weak fields, the relevant degrees of 
freedom are the Goldstone modes of the condensate. These can be studied 
by using a low-energy chiral effective Lagrangian. This leads to a pion-loop 
calculation, which we review in the appendix, that gives a magnetization 



for \eB\ -C mi 



e^B eB 

i^log;;^ for|eS|»m2. 

Using a similar approach and working in the chiral limit, Shushpanov and 
Smilga find that the overall magnitude of the quark condensate is en- 
hanced by a factor 

S(0) ^ 167r2F2 + • • ■ ■ ^'''> 

As the field increases, higher order terms in the chiral Lagrangian become 
important. In any case, the chiral Lagrangian must be abandoned in favor 
of a description in terms of quarks for eB > (47rF^)^ 1 GeV^. It would be 
desirable to understand this transition region between the weak and strong 
field regimes. 

Within the strong field regime, our quark model calculations give an es- 
timate for the magnetization, given in Eqs. (|l^) and i\l7\j . One would like a 



clearer understanding of how the interplay of electromagnetic and QCD ef- 
fects determines the density of quark pairs, thus leading to an improvement 
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upon these estimates. The development of improved techniques for perform- 
ing precise calculations in this regime is clearly needed. Finally, the relation 
between our methods and the more field theoretic approaches to strong field 
chiral symmetry breaking of should be better understood. 
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A Weak-field results 



The response of the QCD vacuum to a weak magnetic field can be com- 
puted by using a chiral effective Lagrangian 0. One can integrate out the 
matter fields to obtain an effective electromagnetic action 5*65 (B). At lead- 
ing order the matter contribution to this effective action, S'^'^**^''(B), arises 
from a single pion loop. Thus, we consider a complex scalar field coupled to 
electromagnetism with action 



S 



d X 



m 



(20) 



Schwinger's classic calculation for a uniform magnetic field then gives 

i\ogdet{-V^V + m^) (21) 



Cmatter 



—I 



ds 



Tre"*^("^''^''+'' 



167r2 



esB 



sinh esB 



(22) 
(23) 



Expanding this in powers of 5, one finds both a quartic and a logarithmic 
divergence. The former is a contribution to the vacuum energy, while the 
latter is the one pion loop contribution to the renormalization of the electric 
charge. Subtracting these divergences gives the renormalized matter effective 
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action 

where e is the renormahzed electric charge. 
The magnetization M = B — H is given by 

A Q X cmatter 

M = B + ^ = ^^^. (25) 
SB SB ^ ^ 

Substituting the result from Eq. ( p^) then leads to 



M = 



dB ] e^B eB r , ^, 2 

— — rlog — - tor eiJ > m 

487r^ m"' 

Since M > 0, the QCD vacuum is paramagnetic. 



(26) 



The effects of higher order terms in the chiral Lagrangian have been 
studied \TE]. These become large for eB ~ (47rF^)^ ^ 1 GeV^, at which 
point the chiral Lagrangian approximation breaks down. For such fields the 
effective photon coupling 



1 _ 1 

4s e2 



q2 n matter 



(27) 



- \- 7^ log for > rn'^ (28) 

is still small, thus justifying our neglect of the quantum fluctuations in the 
electromagnetic field. 
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